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We have obtained the zero-temperature phase diagram of the kagome antiferromagnet with 
Dzyaloshinskii-Moriya interactions in Schwinger-boson mean-field theory. We find quantum phase 
transitions (first or second order) between different topological spin liquids and Neel ordered phases 
(either the V3 x >/3 state or the so-called Q = state). In the regime of small Schwinger-boson 
density, the results bear some resemblances with exact diagonalization results and we briefly dis- 
cuss some issues of the mean- field treatment. We calculate the equal-time structure factor (and 
its angular average to allow for a direct comparison with experiments on powder samples), which 
extends earlier work on the classical kagome to the quantum regime. We also discuss the dynamical 
structure factors of the topological spin liquid and the Neel ordered phase. 

PACS numbers: 75.10. Jm,75.40.Mg,75.50.Ee 



I. INTRODUCTION 

The Dzyaloshinskii-Moriya interactions^ are inevitably 
present in S = 1/2 magnetic oxides when the mag- 
netic bonds lack inversion centers, which is the case 
of the kagome lattice. Although small in strength (it 
originates in the spin-orbit coupling), such a correction 
may favour other phases than the ones usually predicted 
by using the standard Heisenberg model. By break- 
ing explicitly the spin-rotation symmetry of the sys- 
tem, the Dzyaloshinskii-Moriya forces tend to reduce 
the spin fluctuations and may therefore be detrimen- 
tal to the long-searched spin-liquid phases. However, 
would the Heisenberg phase be gapped, such as in va- 
lence bond crystals (generalized spin-Peierls states) or 
in topological spin liquids ^2*2*4 then it would be robust 
against perturbations typically smaller than the gap. 
An example is given by the Shastry- Sutherland com- 
pound SrCu2 (B03)2 which remains in the singlet phase 
in the presence of Dzyaloshinskii-Moriya couplings.— In 
the kagome antiferromagnet (with pure nearest neighbor 
Heisenberg interaction), the very existence of a gap re- 
mains an open question^ algebraic spin liquid and gap- 
less vortex phases have been proposed as alternatives in 
the recent years Current numerical estimates of the 
gap provide a small upper bound ~ 0.05 Jj ^ 10 ! 11 ! 12 ! 13 In 
any case the gap (if it exists) could be smaller than the 
Dzyaloshinskii-Moriya coupling (especially in copper ox- 
ides where it is typically ~ 0.1 J) and the latter is there- 
fore particularly relevant. Experiments on the kagome 
compound ZnCu3(OH)6Cl s 14 T 15 i 16 T 17 have found no spin 
gap (despite a temperature much lower than the upper 
estimation of the gap), but the chemical disorde r 18 ! 19 ! 20 
and, precisely, the existence of Dzyaloshinskii-Moriya 
interaction s 21 ! 22 make the spin gap issue not yet clear. 

In fact, the Dzyaloshinskii-Moriya interactions were ar- 
gued to induce long-range Q = 0, 120° Neel order in the 
kagome antiferromagnet: the algebraic spin liquid theory 



predicts the instability at a critical strength D c = Op 3 , 
while there is a finite quantum critical point at D c ~ 0.1 J 
in exact diagonalization results on samples of size up to 
N = 36*21 Since it is clear that there is no Neel order at 
D strictly zero 1Q i 25 and there is Neel order for D large 
enough^i it is tempting to tackle the problem using the 
Schwinger boson representation for the spin operators^ 
Indeed this approach allows in principle to capture both 
topological spin liquid and Neel ordered phases 27 and 
offers a first framework to describe this quantum phase 
transition. The caveat is that the actual Schwinger-boson 
mean- field solution for the S = 1/2, D = kagome anti- 
ferromagnet is already long-ranged ordered, and it is only 
at smaller values of S (which in this approach is a con- 
tinuous parameter) that a disordered spin-liquid phase is 
stabilized. This result may however be an artefact of the 
mean field approach, and it is possible that fluctuations 
not taken into account at this level do stabilize the disor- 
dered phase for the physical spin- 1/2 system. It is there- 
fore interesting to see what phases the Schwinger-boson 
mean-field theory predicts for the kagome antiferromag- 
net perturbed by Dzyaloshinskii-Moriya interactions. 



In section [II] we present the model and the method. 
In section Uni we discuss the phase diagram together with 
general ground-state properties. In section HVl we illus- 
trate the evolution of observables across the quantum 
phase transition from topological spin liquid to long- 
ranged Neel order: the spinon spectrum, the gap, the 
order-parameter and the condensed fraction of bosons. 
We calculate the equal-time structure factor, its powder 
average and briefly compare both to classical calculation 
and experimental results. In section [V] we present the 
dynamical spin structure factor and its behavior in the 
two phases. We also describe how these behaviors emerge 
from the properties of the spinon spectrum. 
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II. MODEL 

We have considered additional Dzyaloshinskii-Moriya 
interactions (DM) to the standard Heisenberg model on 
the kagome lattice 

H= J2 :-/S,.S, • D,,.iS, ■ S,-)] (1) 

<h3> 

where < z, j > stands for nearest neighbours (each bond 
is counted once) and is a quantum spin operator on 
site z, the Dzyaloshinskii-Moriya field Dij is taken to be 
along the z axis and is staggered from up to down trian- 
gles (See Fig. [1]). For this, we work in a rotated frame 
which allows to eliminate the other components (the S$ 
are to be viewed as rotated operators)* 2 ^ In this case, 
the Dzyaloshinskii-Moriya field favors a vector chirality 
along z and the model $1]) has a global U(l) symmetry 
which can be spontaneously broken or not. 




FIG. 1: The Dzyaloshinskii-Moriya field in the model (spin 
rotated frame -see text) . 

The Schwinger boson representation is written 

S* = \(4 h \)°(Z) n i = a\a l + b\b t = 2S, (2) 

where a \ (resp. b\) creates a boson on site i with spin | 
(resp. |), and a are the Pauli matrices. To fix the length 
of the spin, Sf = (ni/2)(m/2 + 1) = S(S + 1), we need 
to have rti = 2S bosons per site. We define the bond 
creation operator 

Mi = \(e"Mb] (3) 

with Oij = D^/ (2 J). A similar approach has been devel- 
opped by Manuel et al£& on the square lattice. With this 
definition and up to small corrections of order Dfj / J, the 
model takes its standard form>2£ 

H = -2 J A h A ij + NzJS 2 /2, (4) 

<h3> 

where N is the number of lattice sites and z = 4 the 
coordination number. Applied to the vacuum of boson, 
A\j creates a superposition of a singlet and a triplet state 
on the ij bond, i.e. the exact ground-state of a single 



bond of Eq. (|4j). In mean- field theory ^ 26 ! 29 the quartic 
Hamiltonian (|4]) is replaced by a self-consistent quadratic 
Hamiltonian with a bond varying order-parameter Aij = 
(Aij) and the constraint rii = 2S enforced on average 
with Lagrange multipliers A^. Up to a constant the mean 
field Hamiltonian reads: 

H MF = - 2 J A %A l3 + A\ 3 A l3 - | Aj| 2 

<hj> 

- ^A,(n,-2^). (5) 

i 

We now restrict our search to solutions that do not break 
the space group symmetry of the Hamiltonian and hence 
could realize spin-liquid states. There are only four 
classes of such states (called Ansatze in refs i 27 ^ 30 and 
in the following), labelled by their projective symmetry 
group ^ or equivalent ly by fluxes around hexagons and 
rhombus ($ H ex,&Rho) = (0,0), (tt,0), (0,tt), or (7r,7r). 
The flux (j) around a loop (u, • ■ ■ *2n) w ^h an even 
number of links is defined by^i 

Ke* = A 12 (-A* 23 ) . . . i 2n -i,2n(-iL,i)- (6) 

It is a gauge invariant quantity. In each of the four 
Ansatze, all Aij have the same amplitude \Aij\ = A and 
are real in a well chosen gauge. Their signs are repre- 
sented on Fig. [2l These four Ansatze are identical to 
those obtained by Wang and Vishwanath for the kagome 
Heisenberg model: they are fully determined by rota- 
tional and translational invariances of the spin Hamilto- 
nian on the lattice. The (0, 0) Ansatz corresponds to the 
y/3 x \/3 state and the (tt, 0) Ansatz to the Q = state, 
originally found by Sachdev^ while the last two involve 
larger unit cells and may as well be relevant for longer 
range interaction or ring exchange j 3Q i 32 

Using the translation symmetry of the mean-field 
Ansatze, the Hamiltonian is Fourier-transformed, 

Hmf = ^q^q^q + NJzA 2 + (2S + 1)NX (7) 
q 

with 

4>\ = ((«lq) f , • • • , Knq) 1 ", &l-q, • • • , &m-q) (8) 
V N / m X 

(the same for the b operators), m is the number of sites 
in the unit-cell: m = 3 for the first two Ansatze (Af q is 
6 x 6), m = 6 for the second two. A given site is defined by 
a sublattice i G [1, m] and unit-cell x indices. The Hamil- 
tonian is diagonalized using a numerically-constructed^ 3 - 
Bogoliubov transformation 

q = Pj* P q = ( ^ ) (10) 
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FIG. 2: The bond order-parameter Aij of the four symmetric 
Ansatze of model Eq. (pQ). Its modulus is a constant for all 
bonds, an arrow from i to j means Aij > 0. The fluxes 
through the hexagon and the rhombus are: in (a) (0,0), (b) 
(7r,0), (c) (0, 7r) and (d) (7r,7r). The unit-cell is shown by 
dashed lines (twice larger for the last two Ansatze ). 



where U a and V a are m x m matrices. 
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FIG. 3: (Color online) Ground state energies of the 
four Ansatze (blue=(0,0), red=(7r,0), green=(0, 7r), 
magenta=(7r,7r)) vs. = D/2J, for S = 0.025, 0.2, 
1/2 (top to bottom). 



H MF = + NJzA2 + ( 2S + v NX ( n ) 

q/i 

where a; q/x is the dispersion relation of the /i = 1, ...,2m 
spinon modes. Each mode is twice degenerate because of 
the time-reversal symmetry. The ground-state |5) is the 
vacuum of the Bogoliubov bosons. At zero temperature, 
A and A are determined by extremizing the total energy, 
subject to the constraints: 

A=\(A ij )\ (n z )=2S (12) 

(the energy is in fact a saddle point, minimum in A and 
maximum in A^4). 



III. PHASE DIAGRAM 

To obtain the phase diagram, the two self-consistent 
equations ([12]) are implemented numerically for each of 
the four Ansatze of Fig. [2j The energies of the different 
Ansatze are shown in Fig. [3] versus = D/(2J), for three 
values of S: 0.025, 0.2 and 1/2. The corresponding full 
phase diagram of the model is shown in Fig. [4j 

Before discussing the predictions of this model, let us 
remark that in the small S limit all these phases can in 
fact be captured by an analytic perturbative expansion in 
term of flux through closed loops^ At small S, the den- 
sity of bosons is small and the constraints (fT2|) imply that 
Aij must be small compared to A. The mean- field energy 
can then be expanded in terms of gauge- invariant prod- 
ucts of bond order parameters Aij along closed loops. 
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FIG. 4: Phase diagram at zero temperature (spin S,6 — 
D/{2J))\ topological spin liquid (TSL) and Neel long-range 
ordered (LRO) phases characterized by their fluxes through 
hexagons and rhombi. For larger S the region of existence 
of the (0, 0) phase shrinks. Dashed lines are the result of a 
perturbative expansion at small S (see text). 



Following Tchernyshyov et al.^ we have computed the 
expansion up to loops of length 16 in order to calcu- 
late the energy difference between the four Ansatze. The 
results of these calculations give the low S phase bound- 
aries (dashed lines superimposed to the exact results in 
Fig. 2]). It is seen on this example that the so-called flux 



expulsion conjecture^ which predicts that the ground- 
state in non frustrated models has zero flux through any 
closed loop does not apply to frustrated problem, where 
the (0, 7r) and (7r, 0) appear as ground-states in an ex- 
tended range of parameters. 

For the sake of clarity, we will discuss the full phase 
diagram (Fig. [4|) and postpone to the next paragraph the 
illustration on the spinon spectrum of the differences be- 
tween topological spin liquids and Neel ordered phases. 
For 5 = 1/2, there is a direct first-order transition be- 
tween the long-ranged Neel ordered V3 x >/3 and Q = 
phases for a finite Dzyaloshinskii-Moriya coupling. This 
finite range of existence of the x \/3 phase shrinks 
with increasing values of the spin, which is fully compat- 
ible with the classical solution^ 

For low 5 values, and increasing Figs. [3] and [H show 
a sequence of first-order transitions between the (0,0) 
(a/3 x a/3 short-range fluctuations), (0,7r), and (7r, 0) 
(Q = 0, short range fluctuations) spin liquid phases. 
(The (7r, 7r) state is always at higher energy and never 
realized.) The (0,7r) state was argued to be stabilized by 
four-spin interactions up to a large critical SM^ It also ap- 
pears here in a small part of the phase diagram for very 
small 5 but first-order transitions prevent its stability at 
larger 5. 

In the absence of Dzyaloshinskii-Moriya coupling, the 
5 = 1/2 results of this approach are qualitatively not 
consistent with exact diagonalisations which point to a 
non magnetic phase. But in the range of parameters 
around 5 ~ 0.2, the (7r, 0) Schwinger boson mean-field 
results are qualitatively not very far from exact diago- 
nalization results: there are short-ranged Q = corre- 
lations in the Heisenberg cas o 12 i 36 and a second-order 
phase transition to 120° Q = Neel order under the 
effect of Dzyaloshinskii-Moriya coupling^ 

As already mentioned, it may be that a theory be- 
yond mean-field, with a better treatment of the con- 
straint shifts the region 5 ~ 0.2 to the physically accessi- 
ble 5=1/2. Indeed in the Schwinger bosons mean- field 
approach it is well known2& that there are large fluctu- 
ations of the number of bosons. As a consequence, the 
square of the spin operator 

( Sl 2 ) = 5(5+l) + ((nf)-(n,) 2 )/4 (13) 

takes a value 3/2 times larger than 5(5 + 1) (at D = 
0)^i The prefactor is even larger in the presence of 
Dzyaloshinskii-Moriya interaction and amounts to ~ 1.75 
for = 0.25. From the physical spin-1/2 point of view, 
the mean field approximation leads to (unwanted) extra 
fluctuations and, on average, the spin length is larger 
than assumed (because of ([13]) ). The SU(2) symmetry 
of the Heisenberg model can be generalized to Sp(AT) 
(which reduces to SU(2) when J\f = 1) by duplicating 
AT times each pair {a^bi) of boson operators: (c^ a ,^ a ) 
where a = 1 • • • J\f is a "flavour" index^ It can be shown 
that the different boson flavors decouple in the limit of 
large A/", leading to J\f uncorrelated copies of the sin- 
gle flavor problem, for which the exact solution is given 




FIG. 5: Spinon dispersions along Y — K — M — Yoi the first 
Brillouin Zone (see Fig. [6] for the definition) for the (jr, 0) 
Ansatz, S = 0.2. Left: = 0, the system has a gap and is a 
topological spin liquid with short range Q = correlations. 
Right: — 0.25, the lower branch becomes gap less at Y in the 
thermodynamic limit and gives rise to the Goldstone mode of 
the long-range Neel order. 



by the present mean- field treatment of the M = 1 model. 
Thus, it is only in the J\f = oo limit of the model that the 
mean field treatment becomes exact and that the fixed 
"spin length" is recovered. However, the mean- field the- 
ory can describe qualitatively the magnetically ordered 
phases and the deconfined Z2 spin liquid phases of the 
SU(2) models 



IV. SPINON SPECTRUM AND QUANTUM 
PHASE TRANSITION FROM A TOPOLOGICAL 
SPIN LIQUID TO A NEEL ORDERED PHASE 

Spin liquid phases (low S). In the low 5 regime, the 
spinon spectrum of Eq. (fTTj) is gapped everywhere. Fig. [5] 
(left) gives a typical example of such a spectrum for the 
Ansatz (7r,0). The spinon band structure is shown in 
the first Brillouin zone, it has a gap of order 0(1) at 
point T indicative of short range Q = correlations. 
This gap does not go to zero in the thermodynamic 
limit: this phase is a disordered topological spin liq- 
uid. Adding a Dzyaloshinskii-Moriya perturbation has 
two effects on the spectrum. It lifts the degeneracy of 
the lower band corresponding to the symmetry breaking 
of the model from 5/7(2) to U(l) and it reduces the gap 
of the spinon mode that has the chirality opposite to that 
of the Dzyaloshinskii-Moriya field. 

Bose-Einstein condensation. With increasing 5 or 
the gap decreases and above a critical spin 5 C (6 ) ), the 
spinon spectrum shows a finite-size gap, which collapses 
with system size N as 0(1/N). Such a spectrum is shown 
in Fig. [5] (right). In the thermodynamic limit, the bosons 
condense in the soft mode (noted </> qo z oC r with a =|, |). 40 
This gives a macroscopic contribution to the total num- 
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ber of Schwinger bosons: 



-y 



(n t )=N = Y, 



i^ii 2 

s 



(14) 

The condensed fraction x qo in the soft mode is of or- 
der Oil) (or equivalently |Vq ^ | ~ V^V). The transi- 
tion to this Bose-Einstein condensed phase corresponds 
to the development of long-range antiferromagnetic cor- 
relations, as can be seen by computing the static struc- 
ture factor: 



S XX (Q) 



3 

47V 



i,3 



-iQ.(Ri-Rj) 



(d\s?sm) (15) 



where is the position of site i and x is an axis in the 
easy plane perpendicular to Dy. 

The difference in static structure factor between topo- 
logical spin liquid and Neel order phase is illustrated in 
Fig. [6l for the (7r, 0) Ansatz across the Bose-Einstein con- 
densation. In the spin liquid phase, the structure factor 
has broad features located at Q = M e (Fig.[6l left) and at 
equivalent reciprocal points (these are the T points of the 
second Brillouin zone). This structure factor looks very 
similar to exact diagonalization and DMRG resultS i 12 i 36 
The features become sharp in the Bose-Einstein con- 
densed phase (Fig. [6l right), where S xx (M. e ) becomes 
proportional to N: 



47V 



S xx (M e ) 



m 



AF 



C st 



(16) 



where m 2 AF is the order-parameter corresponding to 



0,10 



0,05' 




FIG. 6: (Color online) Static structure factor S xx in the ex- 
tended Brillouin zone for the (7r, 0) state (Q = 0) for a lattice 
size N — 1296. In the spin liquid phase TST^o (left, S — 0.2, 
= 0) there are broad features about M e which become peaks 
with divergent intensity in the ordered phase LRO^^ (right, 
5 = 0.2, = 0.25). 

long-range correlations of the 120° Q = Neel type4i We 
have extracted m AF by fitting the numerical results (up 
to N = 1764) to Eq. ([T6|) with finite-size corrections up 
to order 1/N ^ The extrapolation to the thermodynamic 
limit of rriAF, together with the condensate fraction x qo 
and the gap of the soft spinon are given in Fig. [71 for 
S = 0.2. While near the second-order phase transition, 
the extrapolation of the condensed fraction behaves very 




0.05 



FIG. 7: Gap of the soft mode, order-parameter ttiaf and 
condensed fraction, x qo extrapolated to the thermodynamic 
limit for the (7r, 0) Ansatz as a function of — D/(2J) (S = 
0.2). 
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FIG. 8: (Color online) Static structure factor (powder- 
averaged) across the quantum phase transition from the spin- 
liquid to the Neel phase (S = 0.2 and = (red), 0.05 (green), 
0.1 (magenta) and 0.2 (blue)). Inset: S = 1/2, = 0. 



smoothly and vanishes right at the point where the gap 
opens, the extrapolation of the order-parameter gives a 
small shift. We note that tuaf is very small (a few per- 
cent) in this range of parameters, and, therefore, more 
accurate extrapolations would require using larger sys- 
tem sizes close to the critical point. For strong enough 
Neel order, however, the two order-parameters are clearly 
proportional (ttiaf ^ ^q )- 

We have also calculated the static structure factor for 
powder samples (denoted by S(\ Q|)) by averaging ([T5]) 
over all directions of Q. In the (tt,0) spin-liquid phase 
(Fig. [8] at small D), the overall shape is characteristic 
of short-range correlations of liquids. We can compare 
with the calculation of the diffuse scattering for the clas- 
sical spin-liquid kagome antiferromagnet by Monte Carlo 
simulation^ Here the position of the first broad peak is 
at |M e | instead of |K e | (and the second broad feature is 
at \/7|M e |). This simply reflects the difference of short- 
range correlations of the (7r, 0) Ansatz and the y/3 x y/3 
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classical spin-liquid. In addition, compared with clas- 
sical Monte Carlo simulations, we find no intermediate 
shoulder between the two main broad peaks (except for 
a little hump at >/3|M e |), a point which seems in fact to 
be closer to recent experiments on a spin-liquid deutero- 
nium jarosite^ Moreover, since the response is due to 
quantum fluctuations we expect a rather weak sensitiv- 
ity to the temperature up to temperatures of the order 
of a fraction of J. When D increases, we see the devel- 
opment of the Bragg peaks in the ordered phase, which 
increase as the square of the order parameter when we 
go deep into the ordered phase (Fig. [8] (inset)). In the 
ordered phase we can identify two distinct contributions 
to ([T5]) by using the sum rule 

S **(Q) = ±-Jdu>S xx {Q,u>) (17) 

There are the Bragg peaks (uo = 0) and also the inelastic 
collective modes (which we will detail below) which give 
the additional magnetic background scattering (which 
is the only contribution to scattering in the spin liq- 
uid phase). It is noteworthy that the latter is relatively 
strong for low spin (Fig. [8|) and becomes relatively much 
smaller once the order-parameter is large (inset of Fig. [8]). 
In fact the transfer of spectral weight from the magnetic 
continuum background to the Bragg peaks goes as the 
square of the order parameter. Note also that *S(|Q|) 
does not vanish any more for small Q at D / 0, this 
is because in the presence of the anisotropy the ground 
state is no longer an SU(2) singlet. Although the effect 
is small the measurement at small |Q| in the spin-liquid 
phase could help to figure out what the anisotropy is (or 
give an upper bound when the signal is small, see, e.g. 
Ref . l45l) . In the ordered phase, the finite uniform suscep- 
tibility should give a finite contribution at Q = but we 
recall that this contribution is absent for the U(l) singlet 
ground state of Schwinger boson theory. 

V. DYNAMICAL SPIN STRUCTURE FACTOR 

The Schwinger boson approach allows to calculate the 
dynamical response of the system, which is interesting 
both theoretically and for a direct comparison with ex- 
periments. The inelastic neutron cross-section is propor- 
tional to the spin dynamical structure factor 

/+00 
dte™(S%(t)S° Q (0)) (18) 
-OO 

= 27rJ2\{0\S^\p)\ 2 S(co-co p ) (19) 

P 

where a = x,y,z depending on the polarization of the 
incident neutrons, |0) is the ground-state and 

S Q = ^E e " lQ ' R ' xS - ( 2 °) 

xi 




FIG. 9: (Color online) Dynamical structure factors S xx (c[,uj) 
(left columns) and S zz (^oo) (right columns) for the Q = 
Ansatz in the spin- liquid phase, S — 0.1 (top) or with long- 
range order, S = 0.2 (bottom). The system size is N = 
576. 6 is taken to be 0.25 to emphasize the anisotropy of 
S xx (c[j(jj) and S zz (c[ 1 lo). In the Neel phase, the largest peak 
(in blue) is the quasi-elastic peak, next in intensity are the 
magnon branches (in red) (intensities are cut to see the weaker 
continuum) . 

with R^ x the position of the site zx. We use the Fourier 
transform and the Bogoliubov transformation to express 
Sq in terms of quasiparticle operators (fT0|). At zero 
temperature, since |0) is the vacuum of quasiparticles, 
only creation operators are retained. Given that S^ x is 
quadratic in boson operators, we can only create spinons 
by pair. For example, the following term is present in Sq 

^2 _ ^(Q+q)i^q^^-(Q+q)^q' ( 21 ) 

which applied to the left of the matrix element ([T9]h cre- 
ates a pair of spinons \p) with energy uo p = cj(Q +q ^+a;_ qn 
and wave- vector Q. The intensity of the transition is ob- 
tained by computing the product of matrices, such as in 
([2T]h for each pair of modes. 

We now discuss the general features of the dynamical 
spin structure factor and show the results for the Q = 
Ansatz in Fig. [9j 

In the spin liquid phase (S < S c (0)), all spinons are 
gapped and the two-spinon excitations form a gapped 
continuum, the bottom edge of which is given by the 
minimum of ujq+^i -f cj_ qn over all q and all spinon bands 
(Z,n). S xx (Q,u) and S zz (Q,u) are given in Figs. [9] (top). 
In these figures, we have taken D large enough to see the 
anisotropy of the response (and a small S = 0.1 to be in 
the spin liquid phase). When S increases the lower edge 
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of the continuum shifts to zero and its intensity increases 
continuously. 

For S > S c (6) (Figs. E (bottom), S = 0.2, = 0.25), 
due to the soft spinon mode </> qo £ oC r (See Fig. [5] (right)) the 
system enters a Bose-condensed phase: it has long-range 
correlations and low-energy excitations varying as 1/iV, 
at M e . This spinon mode has a singular contribution 
~ y/N [see Eq. ([H]) ]. As a consequence, the intensities 
have different finite-size scaling, depending on whether 
the pair of excited spinons contains the soft spinon or 
not. We can therefore identify three contributions: 

• Elastic peak. This is the peak with the largest in- 
tensity at Q = M e [shown in blue (online) in Fig. [9] 
(bottom left) and cut in intensity to show the other 
excitations]. It comes from a pair of (identical) 
soft spinons (with wave vector qo = 0) in Eq. (j2T|) . 
and so has energy 0(1/N) and intensity propor- 
tional to |/7 qo i^ l/ qo ^ /V^V| 2 ~ N. By integrating 
over all frequencies, only this (zero- frequency) peak 
contributes to the peak of the (equal-time) struc- 
ture factor, S xx (M e ) (Fig. © (right)). We also note 
that the peak is absent in the zz response (Fig. [9] 
(bottom right)), which is expected because the cor- 
relations are long-ranged in the plane only. 

• Magnon branches. There are three magnon 
branches (shown in red (online) in Figs. [9] (bot- 
tom)), two are "optical" modes, the third one, gap- 
less, is the Goldstone mode of the U(l) symmetry. 
The (almost) flat mode is the weathervane mode 
which is always gapped in the Schwinger boson ap- 
proach, contrary to spin-wave theory^ and irre- 
spective of S.— The magnon here consists of a pair 
of the soft spinon and a spinon of wave-vector Q 
[see Eq. (|2T|) ] so that the magnon dispersion is the 
spinon dispersion coq^^ and the intensities are of 
order \U# n Vq uJ>/N\ 2 ~ 1. 

• Continuum. As for the spin-liquid phase, there 
is a continuum of excitations obtained from con- 
tributions in ([2T]) with two spinons both differ- 
ent from the soft mode. Each peak has intensity 
lU^nV^u/VNl 2 ~ 0(1 /AO and the sum of them 
gives a continuum with finite intensity in the ther- 
modynamic limit, which is absent in lowest-order 
spin- wave theory. 

All these excitations contribute to the sum-rule, 
N((S X ) 2 } = Y. q ${du/2it)S xx (q,u), given the dif- 



ferent density of states. Note that, as expected, 
the Dzyaloshinskii-Moriya interaction introduces an 
anisotropy between the in-plane (left column of Fig. [9]) 
and the out-of-plane (right column of Fig. [9|) dynami- 
cal factors. This anisotropy is visible in the spin-liquid 
as well as in the Neel ordered phase. The effect is more 
spectacular in the latter where the Dzyaloshinskii-Moriya 
interaction strongly suppresses the low-energy response 
in the zz channel around the M e point. 

VI. CONCLUSION 

We have obtained the Schwinger-boson mean-field 
phase-diagram for different values of S. The large S limit 
is in agreement with the semi-classical order by disor- 
der mechanism which selects the y/3 x y/3 stated We 
find that this phase remains stable at small anisotropy in 
a region which becomes broader and broader when the 
spin decreases (and hence quantum mechanical effects 
increase). It is therefore possible in principle to observe 
both ordered phases experimentally and first-order tran- 
sitions between them. However, given the small critical 
strength, the Q = phase is more likely to occur in a real 
compound with large enough S and the kagome potas- 
sium jarosite (S = 5/2) offers such an example . 48 ! 49 

We have identified a region of the phase-diagram (5 ~ 
0.2, Ansatz (tt, 0)) which resembles qualitatively to the 
exact diagonalization results for the S = 1/2 system, 
where the Dzyaloshinskii-Moriya interaction induces a 
quantum phase transition between a topological spin liq- 
uid and the Q = Neel ordered phase. This suggests to 
consider smaller values of the "spin" parameter S as pos- 
sibly relevant within the Schwinger-boson theory, given 
that the treatment on average of the constraint leads to 
enhance (Sf ) compared to S(S + 1). Within this frame- 
work, we could calculate observables such as (i) the cross- 
section of diffuse neutron scattering, with the evolution 
from broad diffuse scattering to Bragg peaks across the 
quantum phase transition (Figs. [6] and [8|) (ii) the neutron 
inelastic cross-section which, in addition to Bragg peaks 
and spin waves, shows a broad continuum in both disor- 
dered and ordered phases (Fig. [9j), absent in the lowest- 
order of spin- wave theory. 
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